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Modeling wealth distribution in growing markets 
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We introduce an auto-regressive model which captures the growing nature of realistic markets. 
In our model agents do not trade with other agents, they interact indirectly only through a market. 
Change of their wealth depends, linearly on how much they invest, and stochastically on how much 
they gain from the noisy market. The average wealth of the market could be fixed or growing. 
We show that in a market where investment capacity of agents differ, average wealth of agents 
generically follow the Pareto-law. In few cases, the individual distribution of wealth of every agent 
could also be obtained exactly. We also show that the underlying dynamics of other well studied 
kinetic models of markets can be mapped to the dynamics of our auto-regressive model. 
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Kinetic models have been drawing substantial atten- 
tion as model markets. In these models]!], markets are 
compared with systems of ideal gases, where agents and 
their wealth are considered analogous to the gas parti- 
cles and their energy respectively. Trading between any 
pair of agent is similar to a collision process where en- 
ergy or wealth is shared between agents. Several models 
of both conserved @, i, i, i, @ and open @, 0, 1, 
economic systems have been studied recently, which dif- 
fer mainly in their exchange rules, namely whether col- 
lision is elastic or inelastic, if a fraction or the whole 
energy of a pair is shared between agents, etc. A mini- 
mal model of a closed market is when a randomly chosen 
pair of particles (agents) collide (trade) elastically such 
that the total energy(wealth) of the pair is shared ran- 
domly between the particles (agents). This wealth con- 
serving dynamics naturally predicts Q a Gibb's distri- 
bution of wealth P{x) ~ exp(—(3x) in equilibrium, which 
has been observed in distribution of income-tax return of 
individuals in several countries Il0( . However, the tails 
of the wealth distribution in ll(J and other economic 



systems follow a power-law distribution, as predicted 
by V. Pareto[T3| (known as Pareto-law). In an attempt 
to get a scale free distribution within the frame-work of 
wealth conserving dynamics, a simple model was pro- 
posed by Chatterjee Chakrabarti and Manna (CCM)0], 
where agents contribute only a fraction of their wealth 
for trading depending on their savings propensities which 
differ among agents. Numerical simulations of this model 
strongly suggest that the the distribution of average 
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with 7 



wealth {wi = (xj) \ fo llow P(w) 

Later, exact results[13( show that the tail of the distribu- 
tion is generically scale- free with 7 = 2. For certain typi- 
cal variations of the model one may get 7 7^ 2. Note that 
in these models, it is only the average wealth which show 
a scale free distribution. Wealth of any individual agent is 
distributed about his average following a Gamma-^fceflij 



distribution. 

Although, kinetic models [l| are successful in describ- 
ing basic features of wealth distribution, they do not cap- 
ture the growing features of most realistic markets. Re- 
cently, growing markets are modeled by pouring an extra 
amount of wealth during each trading which is propor- 
tional to the wealth of one @ or bothg agents partici- 
pating in trading. A power-law distribution for rich was 
observed only in 0] where P(w) ~ w^ 1 ' 7 . To have finite 
average wealth, the tail of this distribution can not be 
scale-free; it must be cut off at some finite w. 

In this article, we introduce a minimal model of grow- 
ing markets and show that this class of models generically 
produce a power-law tail in the wealth distribution, in- 
dependent of the details of the market and the trading 
rules. Both static and growing markets having conserving 
or non-conserving dynamics lead to Pareto-distribution 
of wealth, P(w) w~ J with 7 > 2. Kinetic models are 
just a sub-class satisfying conservation of wealth and 
their dynamics could be mapped to the dynamics of our 
model. This exact mapping suggests that wealth con- 
servation is not necessary for the description of markets. 
It also provides a route to capture the exact distribu- 
tion of the fluctuations of wealth of individual agents. 
Finally, these models, being auto-regressive (AR), bridge 
a connection between kinetic models studied recently in 
econophysics and other AR models [l5| of markets studied 
by economists. 

First, the model. Let us take a system of N indepen- 
dent agents i = 1 ... N, whose wealth at a given time t 
is Xi(t). Each agent i, depending on his investment ca- 
pacity < fii < 1, invests a definite fraction of wealth 
HiXiit) in the market. The market stochastically returns 
a net gain (t). Thus, wealth of agent i at time t is 



Xi(t) = (l-w)xi(t-l)+f(t). 



(1) 
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In this minimal model £(t) is taken as a uncorrelated 
positive stochastic variable with probability distribution 
function (PDF) h(£); it does not depend on {xi}. Thus, 
agents may gain or lose from the market. The auto- 
regressive nature of the model that x(t) depends on x(t — 
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1) is clear from ((TJ). 

First we will calculate the steady state of ([1]). Let us 
define an operator B which moves the variables one step 
backward in time, i.e., x(t — 1) = Bx(t). For convenience 
let us take Xi = 1 — fa, which is similar to the savings 
propensity defined in kinetic models p| . Now, ([TJ can be 
written as, 

*(*) = r^XB m = £ A " c(< ~ n) 

71=0 

t 

n=0 

where we have dropped the index i as agents are inde- 
pendent. In the last step we have used the fact that 
is a uncorrelated random variable and that £(n < 0) = 0. 
Thus, the steady state distribution P(x) which is reached 
as t — > oo is the PDF of the stochastic variable 

OO 

x = A"£(n) (3) 
which is just a weighted sum of {£("■)} with weights{A n }. 

m 

Let x m = A™£(n) be the first m terms of (J3|) and their 

distribution be P m (x). From ([2|) and ([3]) it is clear that 
x m = x(t = to). It implies, first that true steady state 
gets contributions from all orders of A". Secondly, P m (x) 
can be considered as the distribution at t = m. 

Since x rn = A m £(m)+£ m _i, P m {x) satisfies a recursion 
relation, 

Pra{x) P m -l(y)h(^)dy. (4) 

The steady state distribution is then P(x) = 
lim P m (x). Clearly, from (0} one can see that 

m — >oc 

P m (0) = for all m > 0. (5) 

Thus in steady state we must have P(x = 0) = 0. Equa- 
tion (O, being independent of the choice of h(£), can be 
used as generic boundary conditions for l|4"|). Secondly, 
it indicates that the steady state distribution is neither 
Gibb's nor Pareto like, where P(x = 0) is finite. 

To proceed further, we need to be more specific, 
namely we need to know Before considering the 

generic growing markets, we consider few examples of 
static markets where average wealth of the market a = 
(£) is fixed. 

• Normal distribution of £ : The first example is 
when fluctuation of the market is normal, i.e, h(£) 
is a Gaussian distribution denoted by £/(ao,oo) 
with mean «o and standard deviation <7o- In this 
case, the steady state distribution P(x) is Q{a,a) 
where 



It is easy to check that Q(a,a) satisfy Eq. (JT|) in 
steady state; i.e, if PDF of x and £ are G(^q,o-o) 
and Q (/i, a) respectively, then PDF of Ax+£ is same 
as PDF of x. Note that agents in this case can have 
negative wealth even though (x) > 0. The negative 
wealth may be interpreted as debt. 

• Exponential distribution of £: In the next exam- 
ple we take = exp(— £). This case is inter- 
esting, because for A = it gives same steady 
state distribution as that of the CC model||, i.e 
P(x) = exp(— x). For non zero A, we need to solve 
the integral equation J5]). Instead we rewrite it as a 
differential equation (which is possible in this case) , 

^Pm(x) = [P m -l(x) ~ P m (x)} , (7) 

where to > 0, and the boundary conditions are 
given by Eq. ([5]). For to = 0, Po(x) = h(x). In 
terms of G m (s), which is the Laplace transform 
(LT) of P m (x), Eq. ([7]) becomes a difference equa- 
tion 

Gm(s) = — L-G^*), (8) 
whose formal solution is 

m—l 

G m (s) = TJ( 1 + A fe s)- 1 G (s). 

k=0 

Again, let us remind that Go(s) is the LT of 
Pq(x) — h(x). Finally, P(x) is the inverse LT of 

k=l 

which can be written as the following series : 

oo 

P(x) = E C m exp(-x/\ m ) 

rn—1 

oo 

where G" 1 = A m J] (1 - A n ~ m ). (10) 

Although Eq. pH]) is an infinite series, first few 
terms are good enough for numerical evaluation 
of the distribution. Terms up to to = n gives 
P n (x), which can be interpreted either as an ap- 
proximation of true steady state distribution P(x) 
to n th order in A or as the distribution at finite 
time t — n. In Fig. [1] we compare P(x) which is 
obtained numerically with the first four terms of 
([10]) for A = 0.4. Note, that P(x) is a Gamma-like 
distribution similar to what has been obtained in 

Baa. 



• Pareto-law : In our model, the wealth distributions 
of individual agents are not simple and depend 
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FIG. 1: P(x) as a function of x, for h(£) = exp(—£). The 
numerically obtained steady state distribution (for A = 0.4) 
of wealth (symbols) is compared with the first four terms of 
P(x) from Eq. (|10[1 . drawn as solid line. 



on their investment capacities //$. Their averages, 
however, follow a power-law. To prove this let us 
define (xi) = W{. In steady state (x(t)) = (x(t—l)). 
Thus, Eq. U) gives 



<0_ 

Mi 



(11) 



Agents in this model differ in their investment ca- 
pacities. In a system of N agents the average num- 
ber of agents having investment capacity fi is Ng(p) 
where g(jJ>) is the distribution of [i. Thus, we can 
write w(fi) — Distribution of w is then 

PH = 9 (,)A = ®!^M. (12) 

aw w z 

A similar argument was used in [i~3| for deriving 
the wealth distribution of CCM model. Although, 
distribution for the rich (large w) is generically 



P{w) 



' , one can obtain 7 > 2 in typical cases. 



For example, if PDF of /1 is g(/i) = /i a /(a — 1) with 
< a < 1, the asymptotic distribution of (p~2]> re- 



sults P{w) 



r , where 7 = 2 



• Growing markets : The kinetic models of markets 
[J, H, [a| are defined with wealth conserving dy- 
namics, which keeps the total wealth of the sys- 
tem constant. In our model, we can easily in- 
corporate the growth feature of the market ( say, 
stock-markets) by introducing explicit time depen- 
dence in the distribution of £. For example, the 
mean (£) = a(t) may vary in time. The distribu- 
tion of wealth P(x, t) will then depend explicitly 
on t. However, in the adiabatic limit, when a(t) 
varies slowly (such that a(t — 1) w a(t)), we have 
P(x,t — 1) » P(x,t). In this limit, thus, P(x,t) 
is identical to the steady state distribution of the 
time-independent model, where £ has an average 
(0=a(r). 

For demonstration, we take h(£) to be an exponen- 
tial distribution with varying average a(t) = t/T. 



FIG. 2: The figure compares the wealth distribution P(x,t — 
T) of growing markets where average wealth a(t) = t/T, with 
that of the static market with a(t) — 1 (line). As expected, 
For slowly growing market (large T), P(x,t — T) compares 
well with the distribution in static market, suggesting that 
'for adiabatically growing markets, instantaneous distribution 
of wealth is independent of the history'. 



In other words, h(£,i) = exp[— x/a(t)]/a(t). From 
the numerical simulations we calculate the distribu- 
tion P(x, T) at t = T for different values of T. Since 
a(T) = 1, P(x,T) is compared with the steady 
state distribution (flT))) . Figure [21 compares the dis- 
tributions for T = 20 and T — 200, which clearly 
suggests that in the quasi-static limit T — * 00 the 
instantaneous distribution depends only on the in- 
stantaneous distribution of £. 

• Annealed X: Another interesting case is when sav- 
ings propensity of agents change in time. This 
is modeled by taking A as a stochastic variable 
distributed, say uniformly in (0,1). Let h(£) = 
exp(— £). The steady state distribution of wealth 
is then P(x) — r 2 (a;) = xexp(— a;), which can be 
proved as follows. If P(x) = ^2(x) then P(u = 
Xx) = exv(— u)\l(§. Thus, PDF of right hand side 
of Eq. (QJ is r 2 (a;)[l7l which is same as the PDF 
of left hand side. This, completes the proof. 

In rest of the article we discuss kinetic models stud- 
ied in the context of wealth distribution and show that 
the dynamics of these models can be mapped to the AR 
model defined in Eq. ([T]). First let us consider the CCM 
model [6|. The main idea in this model is that the agents, 
labeled by i — 1 . . . N, are considered to have savings 
propensities {A^} distributed as g(X). During trading, 
wealth Xi and Xj of two randomly chosen agents i and j 
changes to x[ and x'j respectively such that 



X^ — XiXi ~\~ 1 T'ij 

x'a = XjXj + (1 - r)T ij: 



(13) 



where Tij = (1 — Xi)xi + (1 — Xj)xj, and r is a stochastic 
variable with PDF U(r), uniform in (0,1). The wealth 
conserving dynamics (|13[) can be interpreted as follows. 
Both agents save a fraction (A,, Aj) of their wealth 
and contribute the rest for trading. The total trading 
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FIG. 3: Comparison of wealth distribution of a tagged agent 
k, having Afc = 0.4 in CCM model (symbol) and the CC model 
(line) where savings propensity is identical (A = 0.4) for every 
agent. Average wealth of CC model is taken as wt = 0.198, 
which is the average wealth of the tagged agent. For both 
CCM and CC model, N = 100. In the inset, we compare the 
distribution of "noise" in these systems. 

wealth Tij is then randomly shared between agents i and 
j, A special case of the model {A^ = A} was studied ear- 
lier by Chakrabarti and Chakraborti (CC)[5J]. Explicitly, 
the dynamics of the model is 

x\ = Xx t + [r(l - X)(xi + Xj)] 

$ = Xx j + [(l-r)(l-X)(x i +x j )], (14) 

where a skew wealth distribution was found for A ^ 0. Al- 
though, the exact steady state wealth distribution P(x) 
of CC model is not known, it can be fitted well to a 
Gamma distribution T n (x) = x™ -1 exp(— x) / (n— 1)! with 
n = (1 + 2A)/(1 — A) [la]. However, later analytical stud- 
ies [19| disagree with the Gamma distribution. In this 
article we refer to it as Gamma-like distribution. 

To see why the dynamics of CCM model is same as 
Eq. (JTJ) we look at the steady state wealth distribution 
of a tagged agent k whose savings propensity is A^ = A. 
Wealth of this agent, which fluctuates in time about the 
mean Wk ~ 1/(1 — A), is distributed as P^ CM (x). Now 
we take a system (namely the CC model) of N agents 
with total wealth WkN and all agents having the same 
savings propensity as that of the tagged agent in CCM 
model, i.e, {\i = A}. Numerically calculated steady state 
distribution of wealth PJp C (x) in this case is found to be 
identical to PJp CM (x). As an example, we take a tagged 
agent in CCM having A = 0.4 and compare PJp (x) and 
P^ CM {x) in Fig. [3] An excellent agreement suggests 
that agents in CCM model are independent, unaware of 
savings propensity of other agents. Since each agent in 
CCM model trade with every other agent, it is not sur- 
prising that this model turns out to be mean-field sys- 
tem of non- interacting agents (it is also observed in [14|]). 
which suggests and supports that one can replace (fT¥j) by 
a single equation 

x\ = Xxi + ij, (15) 

where r\ = rT^ is the noise. Note that Xj satisfies the 
same equation because PDF of r is same as that of 1 — r. 



Replacement of the conserving dynamics (fT4|) by a single 
equation (TPS')) which do not have conservation suggests 
that conservation is not important in these systems. 

To emphasize this point further that 'the wealth con- 
serving dynamics can be replaced by a non-conserving 
one similar to (fT])', we consider other kinetic models. In 
a generic wealth conserving dynamics a pair of agents 
interact as follows, 

X^ = \Xi -\- 1]Xj 

x'j = {l-X)x i + (l-r])x j , (16) 

where both n and A are stochastic variables with PDF 
U(x). We will prove, by mapping wealth conserv- 
ing dynamics (|16[) to a non-conserving one, that the 
steady state distribution of this model is in fact P(x) — 
xexp(— x) = ^2{x) (here (x) = 2). The non-conserving 
dynamics for this model is then 

x = Xx + {, (17) 

where the noise £ = nx and x is the wealth of the other 
agent in the conserving dynamics. Both x and x have the 
same distribution in the steady state. If that distribution 
is T 2 (x), the PDF of f is P(£) = exp(-f)[16]. Thus, the 
dynamics of this model is effectively the same as that 
of the annealed A case of the AR model with exponen- 
tial noise studied earlier in this article, where the steady 
state distribution is P{x) = T2(x). We have done nu- 
merical simulation of the conserving dynamics (| 16[) and 
calculated the distribution of £ = r/x, and the steady 
state distribution of wealth P(x). The resulting distri- 
butions are found (see Fig. [3} to be P(£) = exp(— £) and 
P(x) = a;exp(— x), as expected from the non-conserving 
dynamics. Finally, we take the special case of the model 
with r] — A, which is the kinetic model studied in 0], 
where the steady state distribution is P{x) = exp(— x). 
One may write a non-conserving dynamics in this case as 

x' = X(x + x). (18) 

Again, both x and x have the same distribution in steady 
state. If P{x) — exp(-x), then using [l^, [l7| one 
can show that PDF of right hand side of (TT8"|) is same 
as that of the left hand side. These generic examples 
thus strongly suggest that both conserving and non- 
conserving dynamics approaches the same steady state. 

In conclusion, we introduce a simple model which cap- 
tures the growing feature of realistic markets. Agents 
in these models do not involve in direct trading, they 
interact only through the market. Their net gain de- 
pend on how much they invest and how much they gain 
from the market. The market, naturally noisy, is mod- 
eled by a stochastic variable having a specific distribution 
with fixed or varying mean. In our model, return from 
the market is considered to be independent of individual 
agent's investment (a natural extension would be when 
explicitly depends on x(t)). One of our main re- 
sults is that, the average wealth of agents generically fol- 
low Pareto-distribution. We also argue that, when aver- 
age wealth of the market grows adiabatically, the wealth 
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FIG. 4: Wealth distribution P(x) of a generic kinetic model 
(|16|) (dashed line), obtained from simulations (with N = 100 
and (a;) = 2), is compared with the steady state distribu- 
tion P(x) — iexp(-i) (solid line) of corresponding non- 
conserving model (JXTJ . In the inset we compare PDF of noise 
£ for conserving (dashed line) and non- conserving (solid-line) 
models. For the later one P(£) = exp(— £). 



distribution of agents at any given time depends only 



on the instantaneous market. For static markets, exact 
steady state wealth distribution of agents was calculated 
for a few cases. More importantly, the dynamics of usual 
wealth conserving kinetic models studied in econophysics 
as model markets can be mapped to the dynamics of our 
model which does not have conservation. 

Auto regression is a usual technique for economists, 
for study of financial time series. These new models, be- 
ing auto-regressive in nature, build connections between 
standard auto-regressive models and other kinetic mod- 
els of markets. Kinetic models which are believed to be 
the only model explaining Pareto-law for the tail of the 
wealth distribution is not quite correct. In particular 
both, conservation of wealth during each trading and 
global conservation of wealth are not necessary for ob- 
taining Pareto-distribution. Auto-regression is an alter- 
native which is more generic. 
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